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Q_f Abstract 

■ We consider a nonlinear extension of the generalized network flow model, with the flow 

l/^ . leaving an arc being an increasing concave function of the flow entering it, as proposed by 

Truemper [39j and Shigeno |35| . We give a polynomial time combinatorial algorithm for solving 
corresponding flow maximization problems, finding an e-approximate solution in 0{m{m + 
logn)\og{MU'm/e)) arithmetic operations and value oracle queries, where M and U are upper 
bounds on simple parameters. This also gives a new algorithm for linear generalized flows, an 
^ ■ efficient, purely scaling variant of the Fat-Path algorithm by Goldberg, Plotkin and Tardos [12], 

not using any cycle cancellations. 

We show that this general convex programming model serves as a common framework for 
several market equilibrium problems, including the linear Fisher market model and its various 



m 

>. 

■ extensions. Our result immediately extends these market models to more general settings. We 

^\ , also obtain a combinatorial algorithm for nonsymmetric Arrow-Debreu Nash bargaining, settling 

OO ■ an open question by Vazirani (41j . 

cn 

O ! 1 Introduction 

A classical extension of netvi^ork flows is the generalized network Eow model, with a gain factor 
7e > associated with each arc e so that if a units of flow enter arc e, then j^a units leave it. Since 
^ ■ hrst studied in the sixties by Dantzig [5j and Jeweh [2T] , the problem has found many applications 

^ i including financial analysis, transportation, management sciences, see [21 Chapter 15]. 

In this paper, we consider a nonlinear extension, concave generalized Bows, studied by Truemper 



in 1978, and by Shigeno [33] in 2006. For each arc e we are given a concave, monotone 
increasing function Tg such that if a units enter e then Te{a) units leave it. We give a combinatorial 
algorithm for corresponding flow maximization problems, with running time polynomial in the 
network data and some simple parameters. We also exhibit new applications, showing that it is a 
general framework containing multiple convex programs for market equilibrium settings, for which 
combinatorial algorithms have been developed in the last decade. As an application, we also get 
a combinatorial algorithm for nonsymmetric Arrow-Debreu Nash bargaining, resolving an open 
question by Vazirani [51]. We can also extend existing results to more general settings. 

Generalized flows are linear programs and thus can be solved efficiently by general linear pro- 
gramming techniques, the currently most efficient such algorithm being the interior-point method 
by Kapoor and Vaidya [2^. Combinatorial approaches have been used since the sixties (e.g. 
[211 [27] I38j). yet the first polynomial-time combinatorial algorithms were given only in 1991 by 
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Goldberg, Plotkin and Tardos |12] . This inspired a line of research to develop further polynomial- 
time combinatorial algorithms, e.g. [U HU HH ETJ El [El [IZlllHl [32l [33] ; for a survey on combinatorial 
generalized flow algorithms, see [3l]. Despite the vast literature, no strongly polynomial algorithm 
is known so far. Our algorithm for this special case derives from the Fat-Path algorithm in |12j . 
with the remarkable difference that no cycle cancellations are needed. 

Nonlinear extensions of generalized flows have also been studied, e.g. in [H [3], minimizing a 
separable convex cost function for generalized flows. However, these frameworks do not contain our 
problem, which involves nonlinear convex constraints. 

Concave generalized flows being nonlinear convex programs, they can also be solved by the 
ellipsoid method, yet no practically efficient methods are known for this problem. Hence finding 
a combinatorial algorithm is also a matter of running time efficiency. Shigeno [35j gave the first 
combinatorial algorithm that runs in polynomial time for some restricted classes of functions Fg, 
including piecewise linear. It is also an extension of the Fat- Path algorithm in [12]. In spite of this 
development, it has remained an open problem to find a combinatorial polynomial-time algorithm 
for arbitrary concave increasing gain functions. 

Our result settles this question by allowing arbitrary increasing concave gain functions provided 
via value oracle access. The running time bounds for this general problem are reasonably close to 
the most efficient linear generalized flow algorithms. Concave gain functions extend the applicability 
range of the classical generalized flow model, as they can describe e.g. deminishing marginal utilities. 
From the application point of view, another contribution of the paper is extending generalized flow 
techniques to the domain of market equilibrium computations, where this model turns out to be a 
concise unifying framework. 

The concave optimization problem might have irrational optimal solutions: in general, we give a 
fully polynomial-time approximation scheme, with running time dependent on log(^) for finding an 
e-approximate solution. In the market equilibrium applications we have rational convex programs 
(as in |41j): the existence of a rational optimal solution is guaranteed. We show a general technique 
to transform a sufficiently good approximation delivered by our algorithm to an exact optimal 
solution under certain circumstances. We demonstrate how this technique can be applied on the 
example of nonsymmetric Arrow-Debreu Nash bargaining, where the existence of a combinatorial 
algorithm was open [41] . 

In Section [2l we give the precise definition of the problems considered. Thereby we introduce 
a new, equivalent variant of the problem, called the symmetric formulation, providing a more 
flexible algorithmic framework. Section [3| shows the applications for market equilibrium problems. 
Section [5] explores the background of minimum-cost circulation and generalized flow algorithms, 
and exhibits the main algorithmic ideas. We first present our symmetric generalized flow algorithm 
in Section [5] for the special case of linear gains. Based on this. Section [6] gives the algorithm for 
arbitrary concave gain functions. Section [7] adapts these algorithms for the more standard sink 
formulation of the problems. Section [8] considers the case when the existence of a rational optimal 
solution is guaranteed, and shows how the approximate solution provided by our algorithm can be 
turned to an optimal solution. The final Section [9] discusses possible further directions. 

2 Problem definitions 

We define two closely related variants of the linear and the concave generalized flow problem. Let 
G = {V, E) be a directed graph. Let n = \V\, m = \E\, and for each node i G y, let di be the total 
number of incoming and outgoing arcs incident to i. 

We are given lower and upper arc capacities ^, n : — )• M and gain factors j : E ^ on 
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the arcs, and node demands b : V ^ M. By a pseudoBow we mean a function f : E ^ M with 
^ f ^ u. Given the pseudoflow /, let 



Ci := 



(1) 



In the first variant of the problem, called the sink formulation, there is a distinguished sink 
node t G The objective is to maximize et for pseudoflows satisfying Cj > for all i € V — t. 

This differs from the way the problem is usually defined in the literature with the more restrictive 
Cj = for i G F — t, and assuming ^ = 0, 6 = 0. However, this problem can easily be reduced to 
solving the sink formulation, see e.g. [34j. 

The following extension has been proposed by Truemper [39] and Shigeno [35]. On each arc 
ij G E, we are given lower and upper arc capacities i,u : E M. and a monotone increasing 
concave function Tij : [iij,Uij] — )• M U {— oo}; we are also given node demands b : V ^ M. As for 
generalized flows, a pseudoflow is a function f : E ^ R. with £ < f < u. For a pseudoflow /, let 



In the concave sink formulation, we say that the pseudoflow / is feasible, if ej > for sdl i gV — t 
and et > — oo. The objective is to maximize et for feasible pseudoflows. 

Shigeno [35] defines this problem with Cj = if i G V — t, and 6 = and without explicit 
capacity constraints. She also discusses the version with > 0, and gives a reduction from the 
original version to this one. Whereas capacity constraints can be simulated by the functions T^, 
we impose them explicitly as they will be included in the running time bounds. The formulation 
with ej > seems more natural as it gives a convex optimization problem, which is not the case 
for Ci = 0. 

In the sink formulation, the node t plays a distinguished role. It turns out to be more convenient 
to handle all nodes equally. For this reason, we introduce another, seemingly more general version, 
called the symmetric formulation of both problems. Ideally, we would like to find a pseudoflow 
satisfying Cj > for every i V. The formulation will be a relaxation of this feasibility problem, 
allowing violation of the constraints, penalized by possibly different rates at different nodes. 

For each node i G we are given a penalty factor Mi > and an auxiliary variable > 0. The 
objective is to minimize = X^igy MiKi for a pseudoflow / subject to + > for each i £ V. 

The objective Kf is called the excess discrepancy, kj = means > for each i . These 
conditions might be violated, but we have to pay penalty Mj per unit violation at i. 

The sink version fits into this framework with Mj = oo for i ^ t and Mt = 1. However, it 
can be shown that setting finite, polynomially bounded Mj values, the symmetric version returns 
an optimal (or sufficiently close approximate) solution to the sink version, both for linear and for 
concave gain functions. Besides the sink version, another natural setting is when Mj = 1 for all 
i (zV, that is, maintaining > has the same importance for all nodes. 

While the symmetric formulation could seem more general than the sink version, it can indeed 
be reduced to it. For an instance of the symmetric version with graph G = {V,E), let us add a 
new sink node t with an arc from t to every node i (z V with gain factor 1/Mi. Solving the sink 
version for this extended instance gives an optimal solution to the original problem. The reason 
for introducing the symmetric formulation is its pertinence to our algorithmic purposes. 
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2.1 Complexity model 

The complexity setting will be different for generalized flows and concave generalized flows. For 
generalized flows, we aim to find an optimal solution, while in the concave case, only an approximate 
one. For generalized flows, the gain functions are given explicitly as linear functions, while in the 
concave case, the description of the functions might be infinite. To handle this difficulty, following 
the approach of Hochbaum and Shantikumar |18] . we assume oracle access to the Fjj's: our running 
time estimation will give a bound on the number of necessary oracle calls. Two kinds of oracles 
are needed: (i) value oracle, returning Tij{a) for any a G [lij,Uij\; and (ii) inverse value oracle, 
returning a value /3 with a = Tij{j3) for any a G [rjj(£ij), rjj(ttij)]. 

We assume that both oracles return the exact (possibly irrational) solution, and any oracle 
query is done in 0(1) time. Also, we assume any basic arithmetic operation is performed in 0(1) 
time, regardless to size and representation of the possibly irrational numbers. We expect that our 
results naturally extend to the setting with only approximate oracles and computational capacities 
in a straightforward manner. Notice that in an approximate sense, an inverse value oracle can be 
simulated by a value oracle. 

By an e-approximate solution to the symmetric concave generalized flow problem we mean 
a feasible solution with the excess discrepancy larger than the optimum by at most e. An e- 
approximate solution to the sink version means a solution with the objective value et at most e less 
than the optimum, and the total violation of the inequalities > for i (^V — t is also at most e. 

In both cases, we assume that all Mj values are positive integers, and let M denote their 
maximum. 

For generalized flows, we assume all i, u and b are given as integers and 7 as rational numbers; let 
B be the largest integer used in their descriptions. The running time bound will be 0{m?{rn log i? + 
log M) log n) for the symmetric formulation and 0{'rin?{rn + n log n) log B) for the sink formulation. 
This is the same as the complexity bound of the highest gain augmenting path algorithm by 
Goldfarb, Jin and Orlin [16]. The best current running time bounds are 0(m^'^n^ log -B) using an 
interior point approach by Kapoor an Vaidya [23], and 0(m^n log by Radzik [32], that is an 
enhanced version of [16j . 

For the concave setting, we allow irrational capacities as well; in the complexity estimation, we 
will have U as an upper bound on the absolute values on the 6j's, the capacities iij,Uij and the 
^iji^ij), ^ijiuij) values. For each arc ij, let us define Vij = \rij{£ij)\ whenever Tij{iij) > —00 and 
rij = otherwise. Let 

U = max{max{|5j| : i G V},max{\£ij\, \uij\, \Tij{uij)\,rij : ij G E}} . 

For the sink version, we need to introduce one further complexity parameter U* due to difficulties 
arising iiTij{iij) = —00 for certain arcs. Let U* satisfy U < U*, and that et < U* for any pseudoflow 
(it is easy to see that U* = dtU always satisfies this property). We also require that whenever there 
exists a feasible solution to the problem (that is, > for each i £ V — t and et > —00), there 
exists one with et > —U*. If Tjt{ijt) > —00 for each arc jt G E, then U* = dtU also satisfies this 
property. However, we allow —00 values as in the market applications we also have logarithmic 
gain functions. A bound on U* can be given as in Section [8l 

The main result is as follows: 

Theorem 2.1. For the symmetric formulation of the concave generalized flow problem, there exists 
a combinatorial algorithm that finds an e-approximate solution in 0{m{m + nlogn) \og{MUm/e)). 
For the sink formulation, there exists a combinatorial algorithm that finds an e-approximate solution 
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in 0(m(m + n log n) log(C/*m/e)). In both cases, the running time bound is on the number of 
arithmetic operations and oracle queries. 

The starting point of our investigation is the Fat- Path algorithm [12]. The first main contri- 
bution is the introduction of the symmetric formulation. This is a more flexible framework, and 
thus we will be able to entirely avoid cycle cancellation and use excess transportation phases only. 
Our (linear) generalized flow algorithm is the first generalized flow algorithm that uses a pure scal- 
ing technique, without any cycle cancellation. The key new idea here is the way 'A-positive' and 
'A-negative' nodes are defined, maintaining a 'security reserve' in each node that compensates for 
adjustments when moving from the A-scaling phase to the A/2-phase. 

We extend the linear algorithm to the concave setting using a local linear approximation of 
the gain functions, following Shigeno [35]. This approximation is motivated by the technique of 
Minoux [25] and Hochbaum and Shantikumar [18] for minimum cost flows with separable convex 
objectives. 

3 Applications to market equilibrium and Nash-bargaining prob- 
lems 

Intensive research has been pursued in the last decade to develop polynomial-time combinatorial 
algorithms for certain market equilibrium problems. The starting point is the algorithm for com- 
puting market clearing prices in Fisher's model with linear utilities by Devanur et al. [6], followed 
by a study of several variations and extensions of this model. For a survey, see [261 Chapter 5] or 

In the linear Fisher market model, we are given a set B of buyers and a set G of goods. Buyer 
i has a budget mi, and there is one divisible unit of each good to be sold. For each buyer i £ B 
and good j € G, Utj > is the utility accrued by buyer i for one unit of good j. Let n = \B\ + |G| 
and m be the number of pairs ij with Uij > 0. Let f/max = max{C/jj : i £ B,j £ G} and 
R = maxjmj : i E B}. An equilibrium solution consist of prices pi on the goods and an allocation 
Xij, so that (i) all goods are sold, (ii) all money of the buyers is spent, and (Hi) each buyers i buys 
a best bundle of goods, that is, goods j maximizing Uij/pj. 

The equilibrium solutions for linear Fisher markets were described via a convex program by 
Eisenberg and Gale [8] in 1959; the combinatorial algorithms for this problem and other models 
rely on the KKT-conditions for the corresponding convex programs. Exact optimal solutions can 
be found, since these problems admit rational optimal solutions. 

max mi log Zi 

ieB 

Zi<Yl ^y^y G ^ (EG) 

i€B 

Z,X > 

We show that the Eisenberg- Gale convex program, along with all extensions studied so far, falls 
into the broader class of convex generalized flows. Moreover, in all these extension we may replace 
linear or piecewise linear concave functions by arbitrary concave ones, still solvable approximately 
by our algorithm. 
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For the Eisenberg-Gale program, let us define the graph (y,E) with V = B Li G L) {t}. Let 
ji G E whenever j (z G, i € B, Uij > 0, and set Tji{a) = Uija as a linear gain function. Also, 
let it G E for every i G B with Tit{a) = rriiloga. Finally, set bj = —1 for j € G, and 6j = 
for i ^ B. The above program describes exactly the sink version of this concave generalized flow 
instance with fji = Xij for i €z B, j G G and fu = zi. (To formally fit into the model, we may add 
upper capacities Uji = 1 and uu = Ylji^G^ji without changing the set of feasible solutions.) Hence 
our general algorithm gives an e-approximation for this problem. In Section [HI we show that for 
sufficiently small e we can transform it to an exact optimal solution. 

The flexibility of the concave generalized flow model enables various extensions. For example, 
we can replace each linear function Ujia by an arbitrary concave increasing function, obtaining the 
perfect price discrimination model of Goel and Vazirani pT]- They studied piecewise linear utility 
functions; our model enables arbitrary functions (although a rational optimal solution does not 
necessarily exist anymore). 

In the Arrow-DehTeu Nash bargaining (ADNB) defined by Vazirani |41| . traders arrive to the 
market with initial endowments of goods, giving utility q for player i. They want to redistribute the 
goods to obtain higher utilities using Nash bargaining. The disagreement point is when everyone 
keeps the initial endowment, guaranteeing her q > utility. In an optimal Nash bargaining 
solution we maximize X^jg^ log(zj — Cj) over the constraint set in ()EGp . Unlike for the linear Fisher 
model, equilibrium prices may not exist, corresponding to a disagreement solution. A sophisticated 
two phase algorithm is given in [JTj, first for deciding feasibility, then for finding the equilibrium 
solution. 

The convex program for nonsymmetric ADNB can be obtained from the Eisenberg-Gale program 
by modifying the first set of inequalities to Zi < '^j^q UijXij — Ci. In the formulation as a concave 
generalized flow, this corresponds to modifying the 6j = values for z G to 6j = Cj. Hence this 
problem also fits into our framework. From this general perspective, it does not seem more difficult 
than the linear Fisher model. 

Nonsymmetric Nash-bargaining was defined by Kalai |22J. For ADNB, it corresponds to maxi- 
mizing X^iG-B log(zi — Ci) over the constraint set in (lEGh . for some positive coefficients rrii. The 
algorithm in [51] heavily relies on the assumption mj = 1, and does not extend to this more general 
setting, called nonsymmetric ADNB. Finding a combinatorial algorithm for this latter problem 
was left open in PX]. Another open question in [51] is to devise a combinatorial algorithm for 
(nonsymmetric) ADNB with piecewise linear, concave utility functions. Our result generalizes even 
further, for arbitrary concave utility functions, since the linear functions Uija can be replaced by 
arbitrary concave functions. 

Let G = maxcj. In Section [8] we show how our algorithm can be used to find an exact solution 
to the nonsymmetric ADNB problem in time 0{m{m + n log n)(n log (nC/max^) + logC)). The 
running time bound in [41] for symmetric ADNB {R = 1) is 0{n^ log f/max + n'^ log C). 

Let us also remark that an alternative convex program for the linear Fisher market, given 
by Shmyrev [36j, shows that it also fits into the framework of minimum-cost circulations with 
a separable convex cost function, and thus can be solved by the algorithms of Hochbaum and 
Shantikumar [18] or Karzanov and McCormick p4j. Recently, [42] gave a strongly polynomial 
algorithm for a class of these problems, which includes Fisher's market with linear and spending 
constraint utilities. However, this does not seem to capture perfect price discrimination or ADNB, 
where no alternative formulations analogous to [36] are known. 

As further applications of the concave generalized flow model, we can take single-source multiple- 
sink markets by Jain and Vazirani [20], or concave cost matchings studied by Jain |19] . 

A distinct characteristic of the Eisenberg-Gale program and its extensions is that they are 
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rational convex programs. We may loose this property when changing to general concave spending 
constraint utilities. However, for the case when the existence of a rational solution is guaranteed, 
one would prefer finding an exact optimal solution. Section [8] addresses the question of rationality. 
Theorem 18.11 shows that under certain technical conditions, our approximation algorithm can be 
turned into a polynomial time algorithm for finding an exact optimal solution. We shall verify 
these conditions for nonsymmetric ADNB. 

4 Background and overview 

The minimum-cost circulation proble nfl is fundamental to all problems and algorithms discussed 
in the paper. We give an overview in Section [4.11 We present the two main algorithmic paradigms, 
cycle cancelling and successive shortest paths along with their efficient variants. As already revealed 
by early studies of the problem (e.g. [281 138j ). there is a deep connection between generalized flows 
and classical minimum-cost circulations: the dual structures are quite similar, and the generalized 
flow algorithms stem from the classical algorithms for minimum-cost circulations. In Section 14.21 
we continue with an overview of generalized flow algorithms, exhibiting some important ideas and 
their relation to minimum-cost circulations. We also exhibit here the main ideas of our algorithm 
for the linear case. Section [4.31 considers a different convex extension of minimum-cost circulations, 
when the linear cost function is replaced by a separable convex one. We show how the two main 
paradigms extend to this case, using different approximation strategies of the nonlinear functions. 
Finally in Section 14.41 we consider the concave generalized flow problem, discuss the algorithm by 
Shigeno [35] and its relation to algorithmic ideas of the previous problems. We emphasize some 
difficulties and outline the ideas of our solution. 

4.1 Minimum-cost flows: cycle cancefling and successive shortest paths 

In the minimum-cost circulation problem, given is a directed graph G = {V, E) with lower and 
upper arc capacities ^, n :£"—)■ M U {oo}, costs c :£"—)• R on the arcs and node demands h -.V ^M. 
with Xliey = ^- 



f : E ^ M with i < f < u is called a feasible circulation, if Cj = for all i & V. The objective is 
to minimize f for feasible circulations. 

Linear programming duality provides the following characterization of optimality. For a feasible 
circulation /, let us define the residual graph Gf = {V, Ej) with ij G Ej if ij G E and fij < Uij, or 
if ji € E and iji < fji. The first type of arcs are called forward arcs and are assigned the original 
cost Cij, while the latter arcs are backward arcs assigned cost —cji. For notational convenience, we 
will use fij = —fji on backward arcs. Then / is optimal if and only if Ef contains no negative 
cost cycles. This is further equivalent to the existence of a feasible potential tt : V M with 
T^j — < Cij for all arcs ij G Ef. 

Two main frameworks for minimum-cost flow algorithms are as follows. In the cycle cancelling 
framework (see e.g. [21 Chapter 9.6]), we maintain a feasible circulation in each phase, with strictly 
increasing objective values. If the current solution is not optimal, the above conditions guarantee 
a negative cost cycle in the residual graph; such a cycle can be found efficiently. Sending some flow 
around this cycle decreases the objective and maintains feasibility, providing the next solution. 

^We shall use the term 'circulation' to distinguish from other flow problems in the paper. 
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In the successive shortest path framework (see e.g. [21 Chapter 9.7]), we waive feasibihty by 
allowing e, > or < 0; we call such nodes positive and negative, respectively. However, we 
maintain dual optimality in the sense that the residual graph of the current pseudoflow contains 
no negative cost cycles in any iteration (or equivalently, admits a feasible potential). If there exists 
some positive and negative nodes, we send some flow from a positive node to a negative one using a 
minimum-cost path in the residual graph. This maintains dual optimality, and decreases the total 
ei value of positive nodes. 

For rational input data, both these algorithms are finite, but may take an exponential number 
of steps (and might not even terminate for irrational input data). Nevertheless, using (explicit 
or implicit) scaling techniques, both can be implemented to run in polynomial time, and even in 
strongly polynomial time. 

A strongly polynomial version of the cycle cancellation algorithm is due to Goldberg and Tarjan 
[13j . In each step, a minimum mean cycle is chosen. In dual terms, we relax primal-dual optimality 
conditions to ttj — vTj < Cij +e for ij G Ef, with e being equal to the negative of the minimum mean 
cycle value, decreasing exponentially over time. 

Polynomial implementations of the successive shortest path algorithm can be obtained by ca- 
pacity scaling; the most efficient, strongly polynomial such algorithm is due to Orlin |29j . We 
describe here a basic capacity scaling framework by Edmonds and Karp [7]. Instead of the residual 
graph Ef, we consider the A-residual graph E^ consisting of arcs with residual capacity at least 
A (the residual capacity is Uij — fij on a forward arc ij and fji — iji on a backward arc). The 
algorithm consists of A-scaling phases, with A decreasing by a factor of 2 between two phases. In 
a A-phase, we iteratively send A units of flow from a positive node s with > A to a negative 
node t with et < —A on a minimum-cost path in E^. The A-phase finishes when this is no longer 
possible, which means the total positive excess is at most nA. 

In the A-phase, ttj — tTj < Cij is maintained on arcs of the A-residual graph. When moving 
to the A/2 phase, this might not hold anymore, since the A/2-residual graph contains more arcs, 
namely the ones with residual capacity between A/2 and A. At the beginning of the next phase, 
we saturate all these arcs, thereby increasing the positive excess to at most (2n -|- m)A/2. This 
guarantees that the next phase will consist of at most (2n -|- m) path augmentations. 

4.2 Linear generalized flows — cycle cancefling and excess transportation 

In what follows, we consider the sink version of the generalized flow problem, with sink t gV. For 
a pseudoflow / : E' — >■ R, let us define the residual network Gj = (y,Ef) as for circulations, with 
gain factor jij = l/'y ji on backward arcs. Consider a cycle C in Ef. We can modify / by sending 
some flow around C from some i E V{C). This leaves ej unchanged if j ^ i, and increases by 
(7(C) — 1)q, where 7(C) = I[f,&cle- If 7(C) > 1 then we call C a Row-generating cycle, while for 
7(C) < 1, a Row- absorbing cycle, since we can generate or eliminate excess at an arbitrary node 
i € C, respectively. The amount of flow that can be generated is of course bounded by the capacity 
constraints. 

To augment the excess of the sink t, we have to send the excess generated at a flow-generating 
cycle C to t. Hence we call a pair {C,P) a generalized augmenting path (GAP), ii (a) C is a 
flow-generating cycle, i € V{C), and P is a path in Ef from z to t; or (b) C = 0, and P is a path 
m Ef from some node i with Cj > to t. Clearly, an optimal solution / may admit no GAPs. This 
is indeed an equivalence: / is optimal if and only if no GAP exists. 

The gain factors 7e play a role analogous to the costs Ce for minimum-cost circulations. Indeed, 
C is a flow generating cycle if and only if it is a negative cost cycle for the cost function Cg = 
— log7e. The dual structure for generalized flows is also analogous to potentials. Let us call 
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/i : y — )• M>o U {00} with with /x^ = 1 a label function. Relabeling the pseudoflow / by /i means 
dividing the flow on each arc ij going out from i by /ij. We get a problem equivalent to the original 
by replacing each arc gain by 7^- = jijfii/fij. The labeling is called conservative if 7^- < 1 for all 
ij G Ef, that is, no arc may increase the relabeled flow. 

Assume we have a conservative labeling fi so that = whenever i (z V — t, fii < 00. Let 
V' Q V denote the set of nodes from which there exists a directed path to t. It follows that 
(i) fii < 00 for all i G V, and (ii) V contains no flow-generating cycles. Consequently, given a 
conservative labeling, no GAP can exist, and the converse can also be shown to hold. Note that 
on V' , vTj = — log fii is a feasible potential for Ce = — log7e if and only if fi is conservative. 

Based on this correspondence, minimum-cost circulation algorithms can be directly applied for 
generalized flows as a subroutine for eliminating all flow-generating cycles. This can be indeed 
implemented in strongly polynomial time, see |3H I34|. The novel difficulty for generalized flows 
is how to transport the generated excess from various nodes of the graph to the sink t. In the 
algorithm of Onaga [28], flow is transported iteratively on highest gain augmenting paths, that is, 
from i £ y with ej > on an i — i path P that maximizes 7(i-') = HeePle- It can be shown that 
using such paths does not create any new flow generating cycles. Thus after having eliminated 
all type (a) GAPs, we only have to take care of type (b). Unfortunately, this algorithm may run 
in exponential time (or may not even terminate for irrational inputs). This is due to the analogy 
between Onaga's algorithm and the successive shortest path algorithm - observe that a highest 
gain path is a minimum-cost path for — log 7e . 

The first algorithms to overcome this difficulty and thus establish polynomial running time 
bounds were the two given by Goldberg, Plotkin and Tardos [12] . One of them, Fat-Path, uses 
a method analogous to capacity scaling. A path P in Ef from a node i to t is called A-fat, if 
assuming unlimited excess at i, it is possible to send enough flow along P from z to i so that et 
increases by A. 

The algorithm consists of A-phases, with A decreasing by a factor of 2 for the next phase. In 
the A-phase, we first cancel all fiow generating cycles. Then, from nodes i with Cj > 0, we transport 
flow on highest gain ones among the A-fat paths. This might create new flow-generating cycles to 
be cancelled in the next phase. Nevertheless, it can be shown that at the beginning of a A-phase, 
ej' — et < 2(n + m) A for the optimum value and thus the number of path augmentations in each 
A-phase can be bounded by 2(n + m). Arriving at a sufficiently small value of A, it is possible to 
obtain an optimal solution by a single maximum flow computation. 

The basic framework of [28] and of Fat- Path, namely using different subroutines for eliminating 
flow-generating cycles and for transporting excess to the sink has been adopted by most subsequent 
algorithms, e.g. [14 ^ 116 1 I^T j [Qj [32]. Among them, |16] is an almost purely scaling polynomial time 
algorithm, but it still needs the an initial cycle-cancelling as in [28]. 

In contrast, our algorithm does not need any cycle-cancelling, and adapts Fat-Path to a pure 
successive shortest paths framework. The successive shortest paths algorithms for minimum-cost 
circulations start with an infeasible pseudoflow, having both positive and negative nodes. To use an 
analogous method for generalized flows, we have to give up the standard framework of algorithms 
where > is always maintained for all i € V — t. This is the reason why we use the more flexible 
symmetric model: we start with possibly several nodes having Cj < 0, and our aim is to eliminate 
them. An important property of the algorithm is that we always have to maintain /ij = 1/Mj for 
Cj < 0; for this reason we shall avoid creating new negative nodes. 

Similarly to Fat-Path, we use a scaling algorithm. In the A-phase, we consider the residual 
graph restricted to A-fat arcs, arcs that may participate in a highest gain A-fat-path, and maintain 
a conservative labeling fi with 7^- < 1 on the A-fat arcs. When moving to the A/2-phase, this 
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condition may get violated due to A/2-fat arcs that were not A-fat. Analogously to the Edmonds- 
Karp algorithm, we modify the flow by saturating each violated arc and thereby restitute dual 
feasibility. However, these changes may create new negative nodes and thus violate the condition 
Hi = 1/Mi for ej < we must maintain. 

We resolve this difficulty by maintaing a 'security reserve' of diAfii in each node i {di is the 
number of incident arcs). This gives an upper bound on the total change caused by restoring 
feasibility of incident arcs in all subsequent phases. We call a node A-positive if > djA/ij, 
A-negative if < diAfii and A-neutral if Cj = djA/Xj. A-negative nodes may become negative 
(e-j < 0) at a later phase, and therefore we maintain the stronger condition /ij = 1/Mj for them. 
We send flow from A-positive nodes to A-negative and A-neutral ones. Thereby we treat some 
nodes with Cj > as sinks and increase their excess further; however, as A decreases, such nodes 
may gradually become sources. 

For the sink version, described in Section [TJ we perform this algorithm with Mj = -|- 1 if 
i ^ t and Mt = 1. We shall show that this returns an optimal solution. We remark that the highest 
gain path algorithm [16] can also be modifled to a purely scaling algorithm using the symmetric 
formulation, that enables to start from an arbitrary non- feasible solution and thereby eliminate the 
initial cycle-cancelling phase. 

4.3 Minimum-cost circulations with separable convex costs 

A natural and well-studied nonlinear extension of minimum-cost circulations is replacing each arc 
cost Ce by a convex function Cg. We are given a directed graph G = {V,E) with lower and upper 
arc capacities i,u : E ^ convex cost functions Ce ■ [ie,Ue] — >■ M on the arcs, and node demands 
6 : y — 7- M with Ylii^v ^« ~ ^- ^™ ^° minimize ^^eE ^e(/e) for feasible circulations /. This 
is a widely applicable framework, see [21 Chapter 14]. 

This is a convex optimization problem, and optimality can be described by the KKT conditions. 
Let C^{a) denote the left derivative of Ce- As before, for a feasible circulation / define the 
auxiliary graph Gf = {V,Ef). Let Cij denote the original function if ij is a forward arc and let 
Cji{a) = Cij{—a) on backward arcs. / is optimal if and only if there exists no cycle C in Ej with 
J2eeC ^eife) < 0- ^ual terms, / is optimal if and only if there exists a potential vr : y ^ M 
with TTj — TTi < C^{fij) for all ij G Ef. 

Both the minimum mean cycle cancellation and the capacity scaling algorithms can be naturally 
extended to this problem with polynomial (but not strongly polynomial) running time bounds. 
However, these two approaches relax the optimality conditions in fundamentally different ways. 
The results |24j and jl8] address much more general problems: minimizing convex objectives over 
polyhedra given by matrices with bounded subdeterminants. 

Cycle cancellation was adapted by Karzanov and McCormick [24] . The algorithm subsequently 
cancels cycles in Ef with minimum mean value respect to the C^ife) values. The only difference is 
that the fiow augmentation around such a cycle might be less than what residual capacities would 
enable, in order to maintain 

vr, -vr, <C+(/y) + e ^ijeEf (2) 

for the current potential vr and scaling parameter e. 

For capacity scaling, Hochbaum and Shanthikumar [18] developed the following framework 
based on previous work of Minoux [25] (see also [2i Chapter 14.5]). The algorithm consists of 
A-phases. In the A-phase, each Ce is linearized with granularity A. 

Let Ef' denote the A-residual network. We will maintain A-optimality, that is, there exists a 
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potential tt such that 

^^-^.<92iIll±^tLMM yijeEf. (3) 
A ■' 

Let ^a(u) denote the quantity on the right hand side. If we increase flow on some arc ij by A on 
some ij for which equaUty holds, the resulting pseudoflow will remain A-optimal. We will always 
send A units of flows from a node s with > to a node t with < on a minimum-cost path 
in with respect to 0^{ij). By the above observation, this maintains A-optimality. 

When moving to the next scaling phase replacing A by A/2, we change to a better linear 
approximation of the Ce's. Therefore, ([3]) may get violated not only because contains more 

arcs than E^ does, but also on arcs already included in E^. Yet it turns out that modifying 
each fij value by at most A/2, ([3]) can be re-established. This creates new (positive and negative) 
excesses of total at most mA. 

Recently, |42j gave a strongly polynomial capacity scaling algorithm for a class of objective 
functions, including convex quadratic objectives, and Fisher's market with linear and with spending 
constraint utilities (based on Shmyrev's formulation). Let us also remark that the results [23] and 
[18j actually address much more general problems: minimizing convex objectives over polyhedra 
given by matrices with bounded sub determinants. The framework of [18] needs weaker assumptions 
on the objective function and on the oracle. 

4.4 Concave generalized flows 

As we have seen, both the cycle cancelling and capacity scaling approaches for minimum-cost circu- 
lations naturally extend to separable concave cost functions. Similarly, our algorithm in Section [6] 
for concave gain functions is a natural extension of the generalized flow algorithm in Section [5l 

Nevertheless, we were not able to extend any previous generalized flow algorithm for concave 
gains. Shigeno's [35j approach was to extend the Fat-Path algorithm of [12]. However, [35] obtains 
polynomial running time bounds only for restricted classes of gain functions. The algorithm consists 
of two procedures applied alternately similar to Fat-Path: a cycle cancellation phase to generate 
excess on cycles with positive gains, and a path augmentation phase to transport new excess to 
the sink in chunks of A. For both phases, previous methods naturally extend: cycle cancelling is 
performed analogously to ^24j, whereas path augmentation to [1^. Unfortunately, fitting the two 
different methods together is problematic and does not yield polynomial running time. 

The main reason is that the two approaches rely on fundamentally different kinds of approxima- 
tion of the nonlinear gain functions. While for generalized flows, a cycle cancelling phase completely 
eliminates flow generating cycles, here we can only get an approximate solution allowing some small 
positive gain on cycles at termination. In terms of residual arcs, we terminate with a condition 
analogous to ([2|) for concave cost flows. However, the path augmentation phase needs a lineariza- 
tion of the gain functions analogous to ([3]). Notice that for e = 0, ([2]) implies ([3]) for arbitrary A. 
Yet if some small error e > is allowed, then no general guarantee can be given so that ([3]) hold 
for a certain value A. 

For this reason, our goal was to avoid using the two different frameworks simultaneously. It 
turns out that our scaling based linear generalized flow algorithm (outlined in Section [32]) smoothly 
extends to this general setting. We use the local linearization 0'^{ij) of Tij used by Shigeno, an 
analogue of ([3]). In the A-phase, we consider the graph of A-fat arcs, and maintain 0^{ij) < 1 on 
them. 

When moving from a A-phase to a A/2-phase in the linear algorithm, the only reason for 
infeasibility is due to A/2-fat arcs that were not A-fat. In contrast, feasibility can be violated on 
A-fat arcs as well, as 9^{ij) < 1 < 0^/2{ij) happen due to the finer linear approximation of 
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the gain functions in the A/2-phase. Fortunately, feasibihty can be restored in this case as well, 
by changing the flow on each arc by a small amount. 

5 Linear generalized flow algorithm 

In this section, we investigate the symmetric formulation of the generalized flow problem. In 
describing the optimality conditions, we also allow infinite Mi values to incorporate the sink version. 
However, in the algorithmic parts, we restrict ourselves to finite Mj values. 

We describe optimality conditions in Section 15.11 Notion and results here are well-known in 
the generalized flow literature, thus we do not include references. Section [5. 21 introduces A-fat arcs 
and A-conservative labelings, the feasibility framework in the A-phase. Section 15.31 describes the 
subroutine Tighten-label, an adaptation of Dijkstra's algorithm that finds highest gain A-fat 
paths. The main algorithm is exhibited in Section 15.41 Analysis and running time bounds are 
given in Section 15.51 The final step of the algorithm is deferred to Section 15. 6| where we show that 
when the total relabeled excess is sufficiently small, an optimal solution can be found by a single 
maximum flow computation. 

5.1 Optimality conditions 

Let G = {V, E) be a network with lower and upper capacities i, u and node demands h. In the 
sequel, we assume all lower capacities are 0. Every problem instance of symmetric generalized flows 
can be simply transformed to an equivalent one in this form in the following way. For each arc 
ij G £', increase the node demand hi by lij and decrease hj by —^ijiij- Modify the lower capacity 
of ij to and the upper to Uij — lij. 

For a pseudoflow /, we define the residual network Gf = {V,Ef) as follows. Let ij € Ef, if 
ij € E and fij < Uij or if ji E E and fji > 0. The first type of arcs are called forward arcs, while 
the second type are the backward arcs. For a forward arc ij, let 'jij be the same as in the original 
graph. For a backward arc ji, let 7jj = l/^ij- Also, we define fji = —^ijfij for every backward 
arc ji G Ef. For backward arcs, the capacities are £ji = —JjiUji and uji = 0. By increasing 
(decreasing) fji by a, we mean decreasing (increasing) fij by a/^ij. 

Let P = io . . . ik be a walk in the auxiliary graph Ef. By sending a units of flow along P, we 
mean increasing each fi^^i^^^^ by aIiQ<t<hlitit+i- We assume a is chosen small enough so that no 
capacity gets violated. Note that this decreases e^Q by a, increases Ci^ by aU.Q<t<klitit+i i ^"^^ leaves 
the other Ci values unchanged. 

Let G = i^ . . . ik-iio be a cycle in Ef. Sending a units of flow on G from io modifies only ej^, 
increasing by (7(C) — 1)q for 7(C) = He^cle- G is called a Bow generating cycle if 7(C) > 1. On 
such a cycle for any choice of io S V{G), we can create an excess of (7(C) — l)a by sending a units 
around (assuming that a is sufficiently small so that no capacity constraints are violated). 

The pair (C, P) is called a generalized augmenting path (GAP) in the following cases: 

(a) C is a flow generating cycle, io G V{G), t G V is a node with et < 0, and P is a path in Ef 
from io to t (io = t, P = ^ is possible); 

(b) C = 0, and P is a path between two nodes s and t with > 0, et < 0; 

(c) C = 0, and P is a path between s and t with < 0, et < and 7(P) = Ue^cle > Mg/Mf. 
Lemma 5.1. /// is an optimal solution, then no GAP exists. 
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Proof. In case (a), we can send some a > units of flow around C from ip, and then send the 
generated (7(C) — l)a excess from io to t along P. For sufficiently small positive value of a, this 
is possible without violating the capacity constraints and it decreases the excess discrepancy. In 
case (b), we can decrease the excess discrepancy at t while only decreasing a positive excess at s. 
In case (c), although MgKg increases, MtKt decreases by a larger amount. □ 

The dual description of an optimal solution is in terms of relabelings with a label function 
fj. : V ^ E>o U {00}. For each node i G V, let us rescale the flow on all arcs ij & E hy fii: let 
ftj = fij/lJ-i- We get a problem equivalent to the original one with relabeled gains 7^^- = ^ij^i/^j. 
Accordingly, the relabeled demands, excesses, and capacities are ftf = bi/fii, ef = ei/fii, and 
ufj = Uij/ni. A relabeling is conservative, if for any residual arc ij E Ef, 7^^- < 1, that is, no arc 
may increase the relabeled flow. Furthermore, for each i GV, m > 1/Mi is required and equality 
must hold whenever ej < 0. 

We use the conventions 00 • = and 00/00 = 0. Accordingly, if /ij = 00, we define = ef = 0, 
and 'jj- = for all arcs ji G Ej. If ij € -E/, /ij = 00 and fj,j < 00, then 7^^- = 00. Consequently, if fi 
is conservative, then //j < 00 for any i gV for which there exists a path in Ef from i to any node 
t gV with et < 0. Also, if fi is conservative, there exists no flow generating cycle on the node set 
{i : fii < 00}. This is since for a cycle C, 7(C) = Uij^zcHj = '^ij&clij- 

Theorem 5.2. For a pseudoflow f , the following are equivalent. 

(i) f is an optimal solution to the symmetric generalized flow problem. 

(a) Ef contains no generalized augmenting paths. 

(Hi) There exists a conservative relabeling /i with Cj = whenever l/Mi < /ij < 00. 

Proof. The equivalence of (i) and (Hi) is by linear programming duality, with /ij being the reciprocal 
of the dual variable corresponding to the inequality Cj + Kj > 0. (i) implies (ii) by Lemma [5711 

It is left to show that (ii) implies (Hi). If the excess discrepancy is (that is, ej > for all 
i G V), then /i = 00 is conservative. Otherwise, let N = {t : et < 0}. If Ef contains no directed 
path from i G V to N, then let /ij = 00. For the other nodes i G V, let /ij be the smallest possible 
value of l/(7(P)Mt) for 7(P) = HeeP'Je, where P is a walk in Ef starting from i and ending in a 
node t (z N. By (ii), this is well-defined, since all cycles can be removed from a walk P without 
decreasing 'y{P). 

fi clearly satisfies jijfii/nj < 1. We shall prove /ij > 1/Mj for each i gV and /if = 1/Mj for each 
t G N. If ej > 0, then no such path P may exists as it would give a type (b) GAP. Consequently, 
Ci = 00. If ej < 0, then /ij > 1/Mj as otherwise the optimal P defining /ij would be a type (c) 
GAP. Finally, if t G TV, then fit < l/Aff as the gain of the path P = is defined as 1. □ 

5.2 A-conservative labels 

The residual capacity of arc ij Ef is Uij — fij (for a backward arc ij, this is 7jj/ji). In contrast, 
we define the fatness of ij € Ef by Sf{ij) = ^ij{uij — fij) (on backward arcs, Sf{ij) = fji). The 
fatness expresses the maximum possible fiow increase in j if we saturate ij. This notion enables 
us to identify arcs that can participite in fat paths during the algorithm. In accordance with the 
other variables, the relabeled fatness is defined as Sj(ii) = Sf{ij)/ fXj. 

For the scaling parameter A > 0, we define the following relaxation of conservativity. Let 
H : V M>o U {00} be a label function. Recall that dj is the total number of arcs incident to i. A 
node i e y is called A-negative if < djA, A-iieutrai if ef = djA and A-positive if ef > djA. 
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The A-fat graph E^{A) is the set of residual arcs of relabeled fatness at least A: 

E>}{A) = {ij : ij G Ef : s^^{ij) > A}. 

Arcs in E'^[A) will be called A-fat arcs. The labeling /i is called A-conservative, if 7^- < 1 holds 
for every ij € Ej{A), and /ij > 1/Mj for all i (zV. Further, for every A-negative node i, we require 
Mi = ^/Mi. 

Observe that no nodes with fii = 00 may be present in a A-conservative labeling. If A < A', then 
a A-conservative labeling is also A'-conservative, and 0-conservativity is identical to conservativity. 
Let Ex^{f) = ^jgy max{ef , 0} and Ex^^{f) = Yli&v "^^^{-^i ~ diA,Qi} denote the total relabeled 
excess and total modified relabeled excess for A, respecitively. Note that Ex^{f) < Ex^{f) + 2mA. 

Lemma 5.3. Let f be a pseudoflow with a A-conservative labeling fi. Let < A' < A. Then there 
exists a flow f such that f and fi are A'-conservative and Ex^,{f) < Ex^{f) + 3m(A — A') 

Proof. We shall construct / by modifying / on each arc independently. For A-fat arcs, A- 
conservativity guarantees 7^^- < 1. Consider an arc ij € E^{A') — E^{A), that is, 

A'<7i;-K--/|;)<A 

and assume > 1. Let us set fij = min{uij,fij + ^^-^^-^}. Then ij cannot be A'-fat, since 
either ij ^ Ef, or 7,^.(nf^. - ft^) < A'. 

- ffi = ^i'M - ftj) < A - A'. Also, - /(; < A - A' holds because of 7,^- > 1- We also 
have to consider the possibility that ji is A'-fat for /. In this case, conservativity is guaranteed 
since 7^^. = 1/7^^- < 1. 

To complete the proof of A'-conservativity, we show that / has no A'-negative nodes with 
/ij > 1/Mj. As we have seen, both fij and 'yijfij change by at most A — A'. Consequently for every 
i (z V, the total possible change of the relabeled flow on arcs incident to i is di{A — A'). A node 
is nonnegative for A if > djA and for A' if ef > djA'. Therefore, a nonnegative node cannot 
become A'-negative, proving the claim. 

Further, Ex^,{f) < Ex'^{f) + Yliev ^«(A — A'), and on each arc, at most A — A' units of new 
excess is created. This gives Ex^,{f) < Ex^{f) + 3m(A — A'). □ 

The proof also gives a straightforward algorithm for finding such an /. Let Adjust(A, A') 
denote this subroutine. In particular, Adjust(A,0) finds an / for which /i is a conservative 
labeling. Further, if there are no A-negative nodes for / and fi, then / is a 0-discrepancy optimal 
solution. 

5.3 A-canonical labels 

An edge ij is called tight if 7^- = 1, an a directed path is tight if it consists of tight arcs. Given 
a pseudoflow /, a conservative labeling fi is called canonical, if for each i £ V with /Xj < 00, there 
exists a tight path in Ef from i to a negative node. Analogously, for A > 0, a labeling is called 
A-canonica], if it is A-conservative, and for each i G V there exists a tight path in E^{A) from 
i to some A-negative or A-neutral node. Such a path is a highest gain A-fat path as in [12i . 
(Note that /ij < 00 for every i £ V for A-canonical labelings, and also that paths are allowed to 
end in A-neutral nodes, in contrast to canonical labelings.) By 0-conservative labeling we mean a 
conservative one. 
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Given a A-conservative relabeling fj, which is not canonical, Tighten-Label(/, /i, A) replaces 
^ by a A-canonical labeling /x' with fi'^ > fii for each i &V. 

We first describe Tighten-Label(/, 0), when it is essentially a multiplicative interpretation 
if Dijkstra's algorithm. Let V Q V he the set of nodes i with a directed path in Ef from i to a 
negative node. For nodes in y \ V, let us set /Xj = oo. Let S QV' he the set of nodes i for which 
there exists a (possibly empty) tight path for the current to a negative node. In each step of 
the algorithm, S will be extended by at least one element, and we terminate if S* = V', when the 
current relabeling is canonical. 

If y \ 5 7^ 0, let us multiply fii for each i £ V \ S hy a defined as 

a = min | ^ : ij G ^J/, i G \ 5, j G . 

By the definition of S, a > 1, and after multiplying by a, at least one arc ij G Ej with i G V \S, 
j (z S will become tight. Tight arcs inside S also remain tight, hence S is extended by at most one 
node. Also, the choice of a guarantees that // remains conservative. 

Let us now describe Tighten-Label(/, /i. A) for A > 0. The main difference is that increasing 
fii may turn a A-positive node into A-neutral. We have to stop increasing /.ij at this point and add 
i to S. (This is in accordance with the goal that we allow tight paths to A-neutral nodes as well.) 
In each phase of the algorithm, let S denote the subset of nodes that have a (possibly empty) 
tight path in Ej{A) to a A-negative or A-neutral node. S is initialized as the set of A-negative 
and A-neutral nodes, and is extended by at least one element per phase. The algorithm terminates 
once S = V. 

In every phase, we multiply fii for every z G ^ \ 5" by the same factor a > 1. Consider an arc 
ij G Ef{A) with iGV — S,jGS. This must satisfy < 1. Increasing fii increases 7^^-. Note that 
the fatness s^^{ij) = Sf{ij)/nj is not changed as it is not dependent of f^i. By definition, all nodes 
inV — S are A-positive. When increasing fii, decreases and therefore i may become A-neutral. 
At this point, we have to stop to avoid creating new A-negative nodes. Let us define 



Q = mm < mm 



: U e i?^ ( A) I , min 1 1^ : i G y - 5 [> !> . (4) 



Clearly, a > 1, and after multiplying each /ij by a for i G S", either we obtain a new A-neutral node 
in V — S, or at least one arc ij G Ef with i£V\S,j(zS will become tight. Tight arcs inside 
S also remain tight, and their capacities is unchanged, hence S is extended by at most one node. 
Also, the choice of a guarantees that fi remains A-conservative. Note that arcs with j G V\S may 
disappear from £'^(A) as their fatness decreases. 



5.4 Description of the algorithm 

The algorithm is shown on Figured] We start with /ij = 1/Mj for every node i V , f = and 
A = MB"^ + 1. Once (2n -|- 6m) A < then an optimal solution can be found by a single 

maximum flow computation, as shown in Section [5.61 In this case we terminate. 

The algorithm consists of A-scaling phases. During the A-phase, we maintain a pseudoflow / 
along with a A-conservative labeling /i. The //j values can only increase. Let A^o denote the set 
of A-negative and A-neutral nodes, and D the set of nodes i with > (dj + 1)A. The A-phase 
consists of a sequence of iterations until D becomes empty. In every iteration of the algorithm, we 
update /i to a A-canonical labeling by cahing Tighten-Label(/, //, A). If D 7^ still holds, we 
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Algorithm Symmetric Fat-Path 
for i ^ V do /ij -p-; 
for ij e E do fij ^ 0; 
A ^ MB'^ + 1; 

while (2n + 6m) A > 1/5™ do 
do 

Tighten-Label(/, fi, A); 
D^{i£V -.ay {di + 1)A}; 
No^{ieV :ei< d^A}; 

pick s G D, t £ N connected by a tight path P; 
send A units of flow along P; 
while D / 0; 
Adjust(A, 
A ^ ^• 
Adjust(A,0); 
Tighten-Label(/, /i, 0); 

compute a maximum flow from nodes {s £ V : Cg > 0} 

to nodes {t £ V : et < 0} using tight arcs for fi; 
return optimal primal solution / and optimal dual solution 



Figure 1: The algorithm for symmetric linear generalized flows 



pick an arbitrary s € D, and send A units of flow from s to some A-negative or A-neutral t G A^o 
on a tight path P. At the end of the A-phase, we modify / by Adjust(A, A/2), and proceed to 
the A/2-phase. 

5.5 Analysis 

Claim 5.4. The initial /i is A -conservative, and A-conservativity is maintained during the entire 
A-phase. 

Proof. At the beginning s^{ij) = "fijUij/fij < MB"^ < A for every ij G E, and Ef contains 
no backward arcs. Consequently, E^{A) = and A-conservativity trivially holds. Also, ^ii = 
1/Mi holds for every node i. Tighten-Label(/, /x, A) clearly maintains A-conservativity. This 
is also maintained when sending flow, as we only use tight arcs. At the end of the A-phase, 
Adjust(A, A/2) transforms / to a A/2-conservative pseudoflow. □ 

Claim 5.5. At the beginning of every A-phase, Ex^^{f ) < (2n -|- 3m) A. 

Proof. This holds by definition in the first phase. In the first phase, Ex^{f ) < (X^jgy &«) = 
MB < (2n -|- 3m) A. Once we finish all iterations in the A-phase, D = $ implies Ex^{f) < nA. In 
Adjust(A, A/2), we increase the excess by at most 3mA/2 (Lemma 15.31 ) Hence at the beginning 
of the A/2-phase, Ex'^^^{f) < {2n -\- 3m) A/2, proving the claim. □ 

Lemma 5.6. A A-phase consists of at most 2n + 3m iterations. 
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Proof. Consider the potential function ^ = "^i^y l'<^Six{ef — {di + 1)A,0}/AJ. By Claim [53} 
^ <2n + 3m holds at the beginning. In the relabeling steps, ^ may only decrease, and in every 
path augmentation, it decreases by exactly 1. □ 

Theorem 5.7. The algorithm runs in 0{m{m + n\ogn){m\ogB + log M)) time. 

Proof. A always decreases by a factor of 2, its initial value is MB^ and we terminate if Ex^(f) < 
l/B"^. Hence the number of phases can be bounded by 0(mlogi? + logM). The number of itera- 
tions is 0{m) by Lemma [5.61 The running time of an iteration is dominated by the Tighten-Label 
step, that can be done in 0(m + n log n) time following Fredman and Tarjan's [lOj implementation 
of Dijkstra's algorithm. □ 



5.6 Moving to an optimal solution 

If (2n + 6m)A < l/B"^ at a certain iteration of the algorithm, then by Claim [531 Ex'^{f) + 3mA < 
\/B"^ holds. Next, we transform / to / by Adjust(A,0), so that is a conservative labeling for 
/. By Lemma [Ol Ex'^{f) < l/B^ follows. Then Tighten-Label(/, A) transforms n into a 
canonical labeling. By the following lemma, a single maximum flow computation yields an optimal 
solution. This is the standard technique how most algorithms in the literature terminate. 

Lemma 5.8. Let fi be a canonical labeling for f, and let Gf = {V,Ef) be the subgraph of Gf 
consisting of tight arcs in Ef. If Ex^{f) < then a single maximum flow computation on G 

from source set P = {s € V : eg > 0, < oo} to sink set N = {t (z V : et < 0} terminates with an 
optimal solution. 

Proof. Consider the flow /' resulting after the maximum flow computation. Since flow was sent 
only on tight arcs, fj, is also conservative for /'. If there are no more nodes s with eg > 0, then by 
Theorem 15.21 /' is optimal. Assume now P', the set of such nodes for /' is nonempty. 

Let S* C y be the set of nodes reachable from P' using tight residual arcs in Efi. By optimality, 
S contains no node with negative excess. If an arc ij € E leaves S then either ij is saturated, that 
is, fij = Uij, or 7jj < 1 and flj = 0. Similarly, if ij G E enters S, then = must hold. Also, on 
all arcs ij with i,j G S, /,',• > either 7,^- = 1 or //■ = Uij. Let Fi denote the set of such arcs with 

*J '■J ''J 

fij < "^ij (and 7j^- = 1), and F2 the set of those with ■ = Uij. Therefore, 



ieS ieS \j:ji€E 

= E (rr^f^' - f'/) + E 

Let B* < B^ denote the common denominator of all 7jj's. We claim that every term in the above 
expression is an integer multiple of 1/B*. Indeed, using that fi is a canonical labeling for /, there 
exists a tight path from each node i to a negative node t. l/ni is then the product of the integer 
Mf and the gain factors on such a tight path and is hence an integer multiple of 1/B*. Similarly, 
every "jij/fJ-j is also an integer multiple of 1/B*. Since the Ujj's and 6j's are integers, this verifies 
the claim. Consequently, < Ex^^if) < Ex^^if) < l/B"" < 1/B*. This implies Exf'if) = 0, 
completing the proof. □ 







j:ij€E 




■ E < 


-E^r 


ij£Efn5{S) 
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6 Concave generalized flows algorithm 

We describe the algorithm in the same structure as for the hnear case: Section 16.11 presents the 
optimahty conditions; A-conservative and A-canonical labels are discussed in Sections 16.21 and 16.3^ 
respectively. Section [6.41 presents the algorithm, and Section [6.51 its analysis. 

6.1 Optimality conditions 

The characterization of optimality was given in [35]; we have to modify the results slightly as we 
use the symmetric formulation. Let us call an arc ij € E immense, if Tij{iij) = — oo, and other 
arcs regular. First, let us transform the problem to an equivalent instance with (i) i = for every 
arc and Tij{0) = for every regular; and (ii) every gain function Tij is strictly monotone increasing 
on [0,Uij]. 

For (i), on each arc ij € E, let us replace Uij by Uij — £ij and £ij by 0. If ij is a regular arc, 
we modify the gain function to Tij{a + £ij) — Tij{£ij), and if ij is an immense arc, to Tij{a + £ij). 
Accordingly for every i £ V, let us increase bi by Ylj-ijeE^ij^ ^^"^ decrease it by the sum ofTji{£ji)^s 
on regular arcs. 

For (ii), let us define Uij = mf{p : < p < Uij,Tij ^ijiuij)}- By concavity, T{uij) — 

^ijiuij), and Tij{uij) is strictly monotone increasing on the interval [0,Mij]. Let us replace Uij by 

"^ij ^ Uij . 

For a pseudoflow / : £" ^ R, we define the residual network Gf = (y,Ef) identical as for 
the generalized flow setting: ij G Ef if ij € E and fij < Uij or ji G E and fji > 0. For 
notational convenience, we define fji = —Tij(fij) on backward arcs. We also define the function 

Tji{a) : [-Tij{uij),-Tij{0)] [-Uij,0] by 

r,,(a) = -r-.i(-a). 

Hence Tji{fji) = -fij. 

The concavity of Ty implies that for each < a < Uij, there exists the right derivative, denoted 
by Tfj{a), and for < q < Uij, there exists the left derivative r~.(a). If < A < A', then 

r,,(a + AO-r,,(a) ^ r,,(a + A)-r,,(a) ^ , 

A' - A - 

r,,-(a)-r,,(a-AO ^ r,,(a)-r,,(a-A) ^ 

A' - A - 

for < a < Uij and for A' < a < Uij, respectively. Furthermore, if < a < a' < Uij, then 
rij(a') < r~(Q') < Tfj{a) < T^j{a). The following claim is easy verify. 

Claim 6.1. For any ij G E with < < Uij, T±{fij) = l/TJ^{fji), r-.(^) = l/r±{fj,). □ 

Let P = io . . . ik be a walk in the auxiliary graph Ef. By sending a units of flow along P, we 
mean the following. First we increase fi^i^ by a and set f3 = rjpjj(/jQj^ + a) — rjpj^(/jQjJ to be the 
flow arriving at ii. In step h = 1, . . . ,k — 1, we increase the flow on ihih+i by /3 and set the new 
value of 13 as ^i^ih+iifihih+i + Z^) ~ ^ihih+iifih^h+i)- assume a is chosen small enough so that 
no capacity gets violated. Let /"'^ denote the modifled flow. 

If C = io . . . ik-i is a cycle if Ef, then by sending a units of flow around C from io we mean 
sending a units on the walk io . . . ik-iio- This modifles Ci only in node zq: if the flow increase from 
ik-iio is bigger than a, then Cig increases, and if it is smaller then it decreases. The next lemma 
characterizes when e^Q can increase. For an arbitrary walk P in Ef, let F^ (P) = negpr^(/e). 
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(6a) 
(6b) 



Lemma 6.2. Let C be a cycle in Ef with i G V{C). IfT~j{C) > 1 then Cj can be increased by 
sending some flow around C. IfT^{C) < 1, then it is not possible to increase Cj by sending any 
amount of flow around C . 

Since this property is independent from the choice oii, we simply say that C is a How generating 
cycle if (C) > 1. The lemma is an immediate consequence of the following claim. 

Claim 6.3. Let P = igii ■■■ik be a walk in Ef. For any value of a > 0, the flow increase in ik 
for /"'^ is at most (P)a. On the other hand, for any e > there exists a 6 > so that for any 
< a < 6, /"'^ increases e^^ by at least (F^ (P) — e)a. 

Proof. The first part is trivial by concavity. We prove the second part by induction on the subpaths 
Pfi = iQ . . . ih for h = 0, . . . ,k. There is nothing to prove for h = 0; assume we have already proved 
it for Ph-i- We want to find a 5 > for some an £ > satisfying the claim for Ph. First, it is 
possible to pick a small enough e* > such that 

F+(P,) - e < (F+_,,, J - e*)(F+(P,-i) - e*)- (7) 

By the definition of F^^_^^^, there exists a (5* > such that for any < /3 < 5*, 

By induction for Ph-i and e* , we can choose a small enough 6 > with the following properties: If 
< Q < (5, then F^(P/i_i)a < 6*, and the increase of Ch-i for /"'^^-i is at least /3 = {r^{Ph-i) - 
e*)a. Then ([7]) and ([8]) show that for < a < 5, /"'^ increases e^^ by at least {T~j{P) — e)a. □ 

The definition of GAPs is analogous to the linear case, with the only difference that F^ (P) > 
Mg/Mt instead of 7(P) > Mg/Mt in case (c). The following lemma can be proved similarly as 
Lemma 15.11 

Lemma 6.4. // / is an optimal solution, then no GAP may exist. □ 

Relabelings are also defined analogously as for generalized flows. Given fj, : V ^ M>oU{oo}, let 
us define fj^j = fij/ jii for each arc ij G E. We get problems equivalent to the original with relabeled 
functions F^(a) = Tij{iJLia) / ^j. Accordingly, the relabeled demands, excesses, and capacities are 
6^ = bi/fii, ef = Ci/fii, and ufj = Uij/fii. A relabeling is conservative, if for any residual arc 

ij S Ef, ^ij^ifij) ^ 1) that is, no edge may increase the relabeled flow. Furthermore we require 
> I /Mi for every i and equality whenever < 0. 

We use the same convention for infinite /ij values as for generalized flows. If //j = oo, we define 
bf = ef = 0, ufj = for ij G E, and furthermore Ufi) — ^'^^^ ^ ^f- Finally, for 

ij € Ef with Hi = oo, < oo, let ^f^iftj) = ^■ 

If fj, is conservative, then if for a node i (z V there exists a path from z to a node t € V with 
et < 0, then fii < oo. The following claim is also easy to verify. 

Claim 6.5. F^+(a) = ^T±{a), and F^-(a) = gFr.(a). □ 

This claim implies that for an s — t walk P, Fj^(P) = ^F^ (P), and thus for a cycle C, 
Ff(C7)=F+(C). 

Theorem 6.6 ([35]). For a pseudoflow f, the following are equivalent. 
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(i) f is an optimal solution to the symmetric version. 



(a) Ef contains no generalized augmenting paths. 

(Hi) There exists a conservative labeling fj, with e-j = whenever I /Mi < Hi < oo. 

Proof. The equivalence oi (i) and (Hi) follows by the Karush-Kuhn- Tucker conditions, with /ij being 
the reciprocal of the Lagrange multiplier corresponding to ej + K, > 0. (i) implies (ii) hy Lemma [6.4[ 
The proof of (ii)^(iii) is the same as in Theorem 15.21 with 7(-P) replaced by rjr(P). □ 

6.2 A-conservative labelings 

We define the notion of A-conservative labeling analogously as in Section [5]3] for the linear case.et us 
define the fatness of ij G Ef by Sf{ij) = Tij{uij) — Vij{fij) (if ij is a backward arc, this is equivalent 
to Sf{ij) = fji.) The fatness expresses the maximum possible flow increase in j if we saturate ij. 
This notion enables us to identify arcs that can participite in fat paths during the algorithm. In 
accordance with the other variables, the relabeled fatness is defined as Sj(zj) = Sf{ij)/fj.j. 

Consider a scaling parameter A > 0. The A-fat graph E^{A) is the set of residual arcs of 
relabeled fatness at least A: 

E^f{A) = {ij : ij G Ef : s^(ii) > A}. 

Arcs in E^{A) will be called A-fat arcs. As in [35], we use the following linearization on A-fat arcs 
in chunks of A. 

g^fa-)^= p-Vp , , ^J^E^f{A). (9) 

This is well-defined since Tij(fij) + A/ij < Tij{uij) for A-fat arcs. Q can be written equivalently 
as 

e^i{ij) = — . ij G ^^(A). (10) 

Also, \i ji is a A-fat arc, than using I'jiifji) = — fij and Tj^ = —Tij{—a), we get 

^^^''^ ^ r.,(/^,)-r,^k-A/i.)- ^^^^ 

Consider a label function fi : V ^ M>o U {oo}. A node i G V is called A-negative if ef < djA, 
A-neutral if e-* = diA and A-positive if e-* > diA. The labeling A-conservative, if 9^{ij) < 1 
holds for every ij E E^{A). Furthermore, we require /Xj > 1/Mj for all i (z V, with equality for 
every A-negative node i. 

Using the convexity of r~^, it can be shown that if // is a A-conservative labeling then it is A'- 
conservative for all A' > A. Let Ex^{f) = X^jgy maxjef , 0} and Ex'^{f) = Yliev niaxje^— djA, 0} 
denote the total relabeled excess of positive and A-positive nodes, respecitively. 

The key importance of A-conservativity is that it is maintained when sending A units of flow 
on arcs with 0^{ij) = 1. This is formulated in the next simple lemma. 

Lemma 6.7. Assume fi is A-conservative, and let ij € Ef{A) be an arc with 0^{ij) = 1. If we 
increase f^j by A, then T(fl^j) also increases by A, and A-conservativity is maintained. 
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Proof. 6^{ij) = 1 is equivalent to ^ijiftj) + A) = ^jiftj) + A, showing the first part. Let fij = 
fij + Ani be the modified flow. For the second part, 9^{ij) < 1 for fij easily follows from convexity. 
Further, observe (jlip shows that we get 9^{ji) = 1 for fij. This gives A-conservativity for the 
modified flow as all other arcs are left unchanged. □ 

In contrast to Lemma \b.'6\ the following claim only enables to transform a A-conservative 
labeling can be transformed to a A/2-conservative one. The reason is that besides the set of 
A/2-fat arcs being larger than the A-fat arcs, we may have A-fat arcs with 9^{ij) < 1 < ^^/2(^-^)' 

Lemma 6.8. Let f he a pseudoflow with a A-conservative labeling fi. Then there exists a flow f 
such that 11 is A j 2- conservative for f and Ex'^^^{f) < Ex'^{f) + |mA. 

Proof. Consider a A/2-fat arc ij with 0^i^{ij) > 1 for /, that is, 

r.^' [^ijifij) + - fij < (12) 

There are two possible scenarios: (a) ij was not A-fat, that is, 

yptj < ^ij{uij) - Tijifij) < Afij, (13) 
or (b) ij was also a A-fat arc. Then by A-conservativity, 

r-'(r.,(^) + AM,)-/,, >A^,. (14) 

In both cases, let us define 

fij = (^ijifij) + y/^j) • 

A/2-fatness of ij guarantees that this is well-defined. In case (a), we claim that ij is not A/2-fat 
for /. Indeed, 

Tij{uij) - Tijifij) = Tij{uij) - (^ijifij) + y^j^ < y/ij. 

The last inequality follows by the second part of (jl3p . In case (b), we claim that if ij is a A/2-fat 
arc for / then 0^i^{ij) < 1 must hold for /. Indeed, if we subtract ([12]) from (fH|) . we get 

and by substituting fij, it follows that 

^ij^ (^ijifij) + yMj) - fij > y/^i' 

that is, O^i^iij) < 1 for /. 

We next show that if ji is also a A/2-fat arc for /, then 0^i2{ji) ^ 1 holds for /. Indeed, using 
(fTTI) . O'^^^U^) ^ 1 for / is equivalent to 



^ijifij) ^ij (^fij ~ 2^''' 
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Equivalently, 

A „ A \ A 



^ijifij) + r, ^ij \ fij 2 ^* j — 2^-^' 



Subtracting y/ij, rearranging and applying the monotone increasing function T-^, we get /, 



> 



fij — -jfii, that foUows from (fT2]) by substituting fij. 

We define fij the above way whenever ij is a A/2-fat arc with 9^{ij) > 1. (As a simple 
consequence of concavity, this cannot be the case for both ij and ji.) If this does not hold for 
neither ij nor ji, then let fij = fij. The next claim compares fij and T(fij) to fij and T{fij). 

Claim 6.9. \ft^ - /^.| < # and Irf^-U^.) - rf^{ft^)\ < A. 

Proof. There is nothing to prove if fij = fij. Assume fij was increased as above (the statement is 
equivalent for ij and ji). The first part is identical to (fT2]) . By the definition of fij, 

giving the second part. □ 

For A/2-conservativity, we also need to show that / has no A/2-negative nodes with fii > 1/Mj. 
By the above claim, the total possible change of relabeled flow on arcs incident to i is (i,A/2. A 
node is nonnegative for A if > djA and for A/2 if > (ijA/2. Consequently, a nonnegative 
node cannot become A/2-negative. 

Finally, Ex'^^^{f) < Ex^{f) + ^^^^ (iiA/2, and each arc is responsible for creating at most 
A/2 units of new excess. This gives Ex^^j^if) < Ex^^{f) + yA, as required. □ 

The subroutine Adjust(A) performs the simple modifications described in the proof (in contrast 
to the linear case, this function has only one parameter). 



6.3 A-canonical labelings 

Given a pseudoflow / and a A-conservative labeling the arc ij G E^{A) is called tight if 0^{ij) = 
1. A directed path in E^{A) is called tight if it consists of tight arcs, /x is a A-canonical labeling, 
if from each node i there exists a tight path to a A-negative or to a A-neutral node. Such a path 
is approximately a highest gain A-fat augmenting path. The subroutine Tighten-Label(/, /x. A) 
returns a A-canonical label fJ-' > fJ- for a A-conservative label /i. This is almost identical to the 
algorithm described in Section [5.31 The only difference is in the definition of the multiplier a, which 
is given by ^ for the linear case. Instead, we define 

a = min | min | , . : ij G E't(A) \ , min | — : i G — S* I I . 

In an iteration, we multiply every //j by a for i G V — S, where S is the set of nodes from which 
there exists a tight path to a A-negative or a A-neutral node. We claim that as in the linear case, 
this maintains A-conservativity, and extends S by at least one node. This is a simple consequence 
of the fact that multiplying fii by a multiplies 6^{ij) by a for every incident arc ij. 

To verify that /x remains A-conservative, we also have to check 9^{ij) < 1 on all arcs ij G E^{fi), 
j GV — S. This follows by the convexity of T^-^ . 
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Algorithm SYMMETRIC Concave Fat-Path 
for i G y do /Xj ^ -p-; 
for ij e E do fij ^ mf 
A ^ MU + 1; 
while (2n + 3m) A > e do 
do 

Tighten-Label(/, iJL, A); 
D ^ {ieV ■.€,> {d, + 1)A}; 
iVo ^ G ^ : < d-A); 

pick s & D, t N connected by a tight path P; 

send A units of flow along P; 
while / 0; 
Adjust(A); 
A ^ ^• 

return e-approximate optimal solution /. 



Figure 2: The algorithm for symmetric concave generalized flows 



6.4 The main algorithm 

Let us initialize = 1/Mj for every i & V, and fij = Uij for every ij G E. (We set the upper 
bounds rather than the lower bounds because Tij{0) = —oo is allowed.) Let us pick the initial 
value A = MU + 1. 

The algorithm consists of A-phases, and terminates with an e-approximate solution if (2n -|- 
3m)A < e. During the A-phase, we maintain a pseudoflow / and a A-conservative labeling fi. The 
fii values may only increase. Let D denote the set of nodes i with > {di + 1)A. The A-phase 
consists of iterations, and terminates whenever D becomes empty. In each iteration, we update fi 
to a canonical labeling by calling Tighten-Label(/, /u. A). If Z? 7^ still holds, send A units of 
relabeled flow on a tight path from some s G D to a A-neutral or A-negative node t. 



6.5 Analysis 

Claim 6.10. The initial /i is A-conservative, and A-conservativity is maintained during the entire 
A-phase. □ 

Proof. Initially, f = u and hence Ej is the set of backward arcs. For an arc ij G E, s^{ji) = 
Uij/ni < MU < A, and hence E^{A) = 0. Also, = 1/Mj holds for every node i. Tighten- 
Label(/, /i. A) clearly maintains A-conservativity. We use only tight arcs to send flow, and 
Lemma [621 guarantees that this preserves A-conservativity. At the end of the A-phase, Adjust(A) 
transforms / to a A/2-conservative pseudoflow. □ 

Claim 6.11. The A-phase starts with Ex^{f) < (2n-|-3m-)A. 

Proof. For the initial solution, Ex^{f) < M{'^-^y \bi\ -\- mU) < {m + n)MU. The claim follows, 
since A = MU -\- 1. Once we flnish all iterations in the A-phase, D = % implies Ex^^{f) < nA. By 
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Lemma 16.81 Adjust(A) transfroms / to a A/2-conservative solution by increasing the excess by 
at most |mA. Hence the A/2 phase starts with Ex^^{f) < {2n + 3m)A/2, proving the claim. □ 

Lemma 6.12. A A-phase consists of at most 2n + 3m iterations. 

Proof. Consider the potential function = ^■gy[max{e^ — {di + 1)A,0}/AJ. By Claim [OTI 
^ < 2n + 3m holds at the beginning. In the relabeling steps, ^ may only decrease, and in every 
path augmentation, it decreases by exactly 1. □ 

Recall that kj = J2ieV ■^i'^i ~ X^ieV ^» ™™{~^*' ^1 denotes the excess discrepancy. For 
a A-conservative fi, MiHi = = holds for every node i with Cj < 0, because of //j = 1/Mj. 
Consequently, is the total relabeled deficiency of the negative nodes. The next theorem shows 
that if A < e/{2n + 3m), then we have an e-optimal solution at the end of the A-phase. 

Theorem 6.13. At the end of phase A, the actual f is {2n + 3m) A- optimal. 

Proof. Let us keep running the algorithm forever unless it finds a 0-discrepancy solution at some 
phase. First, consider the case when for some A' = A/ 2^, we terminate with a 0-discrepancy 
solution. In all phases between A and A' , the total decrease of excess discrepancy is bounded by 
{2n + 3m){A/2 + A/A + ...+A/2^) < (2n + 3m)A. Since in the A'-phase we have a 0-discrepancy 
solution, the total discrepancy at the end of the A-phase is at most (2n -|- 3m)A, proving the 
theorem. 

Assume now the procedure runs forever. For each i ^ V , m is decreasing and thus converges 
to some limit n*. Let k* = X^jgyMjK*. As above, the total decrease of the excess discrepancy 
after phase A is bounded by (2n + 3m) A, hence kj < k* + (2n -|- 3m) A. The proof finishes by 
constructing an optimal pseudoflow /* with discrepancy k* . 

Let denote the flow at time t, for A*^*) = Ao/2*, with labels fif\ For each node i, fif^ is 

increasing; let ^* = limt^oo fJ-f^ ■ For every ij G £', f^j^ is a bounded sequence (0 < /^^^ < Uij). 

Consequently, we can choose an infinite set T' C N so that restricted to t € T', all sequences f-j^ 
converge; let f*j denote the limits. We shall prove that /* is an optimal pseudoflow with optimal 
labeling completing the proof. 

Let Voo = {i : fi* = oo}. We claim that V — Voo 7^ 0. Indeed, if i is A-negative in a certain 
phase, then /ij = 1/Mj, and once i becomes A-positive or neutral, it would never again become 
A-negative. Consequently, the set of A-negative nodes is decreasing. Once it becomes empty, we 
arrive at a 0-discrepancy solution. If it never becomes empty, then we have a set A^* which remains 
the set of A-negative nodes after a finite number of steps and thus //^ = 1/Mj for i € N* . 

Let e* denote the excesses of /* If e* < 0, then clearly i G N* and /x* = 1/Mj. If e* > 0, 
we shall prove fj,* = 00. For a contradiction, assume /x* < cxd. Then for sufficiently large t £ T', 
{di + 2n + 3m)A(*Vf^ < ef^ and thus Ex^^^,^{f) > {2n + 3m)A(*), a contradiction. 

We have to prove rf^*+(j^^*) < 1 whenever ij € Ef*. If /i* < 00, then Tf.^{f*f) = 0. If 
fi* < 00, then the definition ([9j) gives 

- ""^'^ r-(r.,(4)) + AW;.f)-4) 

Then AW//j*^ and hence the first fraction converges to ^tjiftj) = '^/{^i/'}^i^tj{fij)), while 
the second to /U* //i* , leading to the conclusion using Claim 16.51 □ 
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Theorem 6.14. The above algorithm finds an e-approximate solution to the symmetric concave 
generalized flow problem in 0{m{m + nlogn)log{MUm/e)) oracle calls. 



Proof. The initial value of A is MU + 1, and we terminate if A < e/{2n + 3m) by Theorem 16.131 
Hence the total number of scaling phases is 0{log{MUm/e)). The number of iterations in a phase 
is 0{m) by Lemma 16.121 and the running time of an iteration is dominated by Tighten-Label, a 
slightly modified version of Dijkstra's algorithm that can be implemented in 0{m + nlogn) time 
using Fibonacci heaps as in [10] □ 



7 Sink versions of the problems 

In this section, we show how the algorithms in Sections [5] and [6] can be applied to solve the to 
the sink versions of the corresponding problems. For generalized flows, let us set Mf = 1 and 
Mi = + 1 for every ieV-t. Let us set bt = [Ej-jieE IjtUjt - Ej:tj€E ^tj + 1] < dtB^ + 1. This 
a strict upper bound on Yl,j:jt<^Eljtfjt — '^j-.tjeE ftj-, hence < will hold for any pseudoflow. 

Let us run the algorithm for the symmetric formulation with these Mj's, returning an optimal 
pseudoflow / and optimal labels /i. We claim that / is also optimal for the sink formulation. If 
Cj > for all i ^ t, this is clearly the case. 

On the other hand, we claim that if there exists a node i with < 0, then the sink version 
is infeasible. First, we show that such an i cannot be reached from i on a path in Ef. Indeed, 
if P were a t—i path in Ef, then 1 > 7^(-P) = j{P)fit/ f^i- Since both i and t are negative, 
Hi = I /Mi = 1/(5" + 1) and ^it = 1. Consequently, 7(P) < 1/(5" + 1). This is a contradiction 
since ^{P) is the product of at most n rational numbers, each with denominator at most B. Let 
y C y be the set of nodes j for which there exists a path in Ef from j to some i — t with 
ej < 0. This set verifies that the sink version cannot be feasible. 

By setting the bt value and the Mj's, B has increased to dtB"^ + 1 and M = i?" + 1. This gives 
running time 0{m?{rn + nlogn) log 5). 

Let us turn to concave generalized flows. An e-approximate solution to the sink version means 
a pseudoflow / with ^^gy_^max{0, — ej} < e and et being at least the optimum value minus e. 

Let us set bt = U* + 1, a, strict upper bound on J2j-jt£E^jt(fjt) ~ J2j-tj€E ftj defined U* 
in Section [2?T]) . Thus ej < is always guaranteed. Let us set Mj = [2[/*/e] + liiiGV — t and 
Mt = 1. Let us run the algorithm for the symmetric formulation to obtain an e-optimal solution /. 

If Kf > 2U* + e, then no feasible solution may exist. Indeed, by the definition of U*, if there 
is a feasible solution /', then there exists one with et > —U* . If /' is such a feasible solution 
for the sink formulation, then its excess discrepancy for the symmetric formulation is at most 
Kf < bi + U* < 2U* , a contradiction as / was e-optimal for the symmetric formulation. 

If Kf < 2U* + e, then 

Also Kt cannot be further than e from the optimum value of et for the sink formulation. Indeed, let 
/' be the optimal solution to the sink formulation with e't flow reaching the sink. Then Kf/ = bt — e't. 
The claim follows by 

bt - e't + e = Kf> + e > Kf > Kt = bt - et, 
and thus et > e[ — e. This gives a running time bound 0{m{m + nlogn) log{U*m/e)). 
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8 Finding the optimal solution for rational convex programs 



In this section, we first give a general theorem which shows how an approximate solution to the sink 
version can be converted to an exact optimal solution, given that one exists. We shall verify the 
required technical properties with appropriate parameters for nonsymmetric Arrow-Debreu Nash 
bargaining. Unlike the linear Fisher model, ADNB might be infeasible. However, it can be shown 
that if the problem is infeasible, then for appropriate (polynomially small) e, the e-approximate 
version is also infeasible. Similar reductions should hold for all other rational convex programs 
discussed in Section [3] as well, giving polynomial time algorithms for finding optimal solutions. 

Theorem 8.1. Let problem V be given by the sink formulation with n nodes and m arcs, and 
complexity parameters U, U* . Assume V is guaranteed to have a rational optimal solution, and the 
following conditions hold for some values e, T and a function T{n,m,U*). 

(PI) Consider the algorithm for the sink version for an e- approximation. Then either there is 
no e-feasible solution, or /Xj < T holds for any i ^ V , even if running the algorithm for an 
arbitrary number of phases. 

(P2) A subroutine is provided for finding an optimal solution f in T{n, m, U*) time, if the following 
assumptions hold. Assume that for each ij € E, we are given an interval lij C [£ij,Uij] with 
< 2T£, with the guarantee that there exists an optimal solution f* with f*- G Iij for all 
ij G E. 

Then there exist an algorithm for finding the exact optimal solution or proving that the problem is 
infeasible in 0{m{m + nlogn) log(C/*m/e)) + r(n, m, U*). 

We remark that in (P2)., / = /* is not required. 

Proof. Let us formulate the symmetric version for e-approximation as in Section [71 Assume we run 
the algorithm for this problem forever, as in the proof of Theorem 16. 131 The /Xj's shall converge to 
some finite values /x* < T. In any A-phase, the total change of f^ - is bounded by e' = (2n + 3m)A, 
and thus fij may change by at most Te' . Therefore all /j^'s converge to some values flp which can 
be seen to give an optimal solution, as in the proof of Theorem 16.131 

The algorithm terminates whenever A < e/(2n + 3m). At this point, the intervals lij = 
[fij — Te,fij+Te] satisfy the conditions in (PI), since \fij — fij\ < Te. Running the e-approximation 
algorithm and then the algorithm in (P2) gives the running time bound. □ 

To ensure property (PI), a useful method is to enforce the existence of a unique optimal 
solution by perturbing the input data, as done by Orlin [30j for linear Fisher markets. If there is 
a unique rational optimal solution /* with all entries having denominator at most Q, then setting 
2Te < 1/Q enables us to identify the set of arcs with f*j > 0. This can be already enough to 
compute /* efficiently. 

8.1 Application to nonsymmetric Arrow-Debreu Nash bargaining 

Let us now apply Theorem 18.11 for the nonsymmetric ADNB problem. Let us assume all utilities 
Uij, budgets m, and disagreement utilities q are nonnegative integers, with C/max = max{f7jj : i G 
B,j G G}, R = max{mj : i G B}, and C = maxjcj : i G B}. Let n = \G\ + \B\ and let m be the 
number of pairs ij with Uij > 0; in the concave generalized flow instance, the number of nodes is 
n + 1 and the number of arcs is m + \B\. Let us assume that there exists at least one arc with 
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positive utility incident to any buyer and to any good. The special case c = is identical to Fisher's 
market with linear utilities. 

Consider a candidate solution with price pj for each good j £ G. Let Xij > denote the amount 
of good j purchased by buyer i. It follows from the KKT-conditions (see also ^T]) that {p,x) is 
an optimal solution if and only (i) YlieB -"^v ~ ^ ^'-'^ each good j, that is, each good is fully sold; 
and (ii) for any buyer i and good j, Uij/pj < {J2jeG ^ij^ij ~ Ci)/mi, and equality holds whenever 
Xij > 0. 

By strict concavity of the objective, the utilities X^jgG ^ij^ij accrued by the players are the 
same in any optimal solution whenever there is a feasible solution. Yet these same values can be 
obtained by different allocations. As in [30], we assume that there exist a unique optimal allocation 
as well. This can be done by a lexicographic perturbation of the Uij values, without significantly 
increasing the running time. This guarantees that the set of arcs with Xij > in an optimal solution 
is cycle-free. 

For the concave generalized flow instance, let us set upper capacity Vji = 2 on each arc ji with 
j G G, i € B, and vu = 2 J2jeG each arc it (we set larger capacities so that the capacity 

constraints would never become tight). Hence the complexity parameter U is bounded by 2|G|i7max- 
We shall prove the following. 

Theorem 8.2. Let K = nRU^ni,^. Setting T = U* = max{C,nKlogK}, e = 1/(2K"C/*) satisfies 
the requirements on U* in Section \2.1\ and (PI) and (P2) in Theorem \8.1\ Our algorithm delivers 
an optimal solution in running time 0(m(m + n log n)(nlog(ni?C/niax) +logC)) for nonsymmetric 
ADNB. 

The running time T{n,m,U*) will be negligible compared to the main algorithm and therefore 
it does not affect the complexity bound. In the rest of the section, we shall verify the above choices. 

Lemma 8.3 (see |4H Thm 2], |30l Lemma 2.1]). Assuming that the problem is feasible and there 
exists a unique optimal allocation x* , all positive x*j values are rational numbers with a common 
denominator S < K^. 

Proof. The optimal allocations x* and prices p* can be uniquely obtained given the set F of arcs 
ji with x*j > 0. If we introduce the variable Qj = 1/pj, then an optimal solution must satisfy the 
following system of linear equations. 

Xij = 1 Vi G G (15) 

j:jieE 

UikXik - UijmiQj = Ci \/ji € F 

k:kieE 

We claim that this system has a unique solution {x* ,p*). Fixing the price of one good jo as pj^ = a 
in a component of F, it uniquely determines pj for any good j in the same component as a times 
the product of the Uab values on the unique path from jo to j. Similarly, this determines the best 
bang-per-buck values bi = Uij/pj = (Ylk-ki&E^ikXik — Ci)/mi for ji G F, which are proportional 
to 1/a. The optimality conditions imply that in an equilibrium, the money spent by buyer i is 
ri = rrii + Ci/bi. In each component of F, the sum of prices should be equal to the money spent by 
the buyers. This uniquely determines all prices and bang-per-buck values in the component. The 
Xij values in the component have to sum up to 1 for each good i and J2k-ki&E ^ikXik = hmi + q. 
As F is a forest, this has a unique solution. 

In the solution to ()15p . a common denominator is the determinant S oi a largest non-singular 
submatrix of the constraint matrix. The Hadamard-bound gives S < (nC/max-R)"' = K"^, using that 
IFI < n - 1. □ 
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The above proof also gives a simple linear time algorithm for finding the optimal solution, 
verifying (P2) if 2Te < l/K"^, with r(n,m,?7*) being negligible compared to the running time of 
the approximation algorithm. 

Next we justify the choice of [/*. U < U* clearly holds. For an arbitrary pseudoflow /, 
et < ^jgQ m-j log tijf < nRlogU < U* . It is left to show that if the problem is feasible, there 
exists a feasible solution with et > —U*. Since the Uij and Cj values are integers, whenever 
^j-ji£E ^ij^ij — Ci>0, it should be at least 1/K". Consequently, if the sink version of the problem 
is feasible, the optimal objective value is at least et > ^^g^ mj log(l/iir"') > —n^i? log > —U*. 

The next claim verifies (PI) and thus completes the proof of Theorem 18.21 

Lemma 8.4. Either the problem is not feasible, or < U* holds for any k £ B U G in arbitrary 
A-phase. 

Proof. Recall from Section[7]that we solve the sink version by reducing it to the symmetric algorithm 
with Mt = 1 and Mi = [2f/*/e] + 1. Since /ij is non-decreasing, these values converge to some 
limits //* E M U {oo}. We have < fij < Uij on all arcs ij £ E in every phase, and therefore we 
can choose an infinite subset T' C N so that all /jj's converge if we restrict ourselves to iteration 
numbers in T'. As in Theorem 16. 13| it can be easily verified that the limit /* is an optimal solution 
to the symmetric version with conservative labeling /i*. 

As in Section [71 if > 2U* , then the sink version is not feasible, and otherwise /* is also 
optimal to the sink version. In the latter case, fl^ > for arbitrary i £ B must hold as otherwise 
ej = — oo in the sink version gives infeasibility. Recall also that the symmetric version was defined 
such that < holds for every feasible solution and therefore fj,* = 1. 

Consider now an arbitrary i £ B. Both it, ti G Ef (it is easy to verify that fa = vu is impossible), 

it* 

and therefore ■j^j^ = l must hold by the conservativity of /i*. This means //* = f*t/mi < U < U*. 

Finally, let j S G. Then for an arbitrary arc ji E E, ji € Ej easily follows, and therefore 
conservativity gives j^Uji < 1, which implies fi* < U*. □ 

Finally, we remark that if we apply this algorithm to linear Fisher markets (c = 0) , the algorithm 
runs in a fundamentally different way as [6j or |30j. While both these algorithms increase the prices, 
our algorithm works the other way around: it starts with the highest possible prices, and decreases 
them. 

9 Discussion 

We have given the first polynomial time combinatorial algorithms for both the symmetric and the 
sink formulation of the concave generalized flow problem. Our algorithm is not strongly polynomial. 
In fact, no such algorithm is known already for the linear case: it is a fundamental open question 
to find a strongly polynomial algorithm for linear generalized flows. If resolved, a natural question 
could be to devise a strongly polynomial algorithm for some class of convex generalized fiow prob- 
lems, analogously to the recent result [42], desirably including the market and Nash bargaining 
applications. 

Linear Fisher market is also captured by [32]. A natural question is if there is any direct 
connection between our model and the convex minimum cost flow model studied in [12] • Despite 
certain similarities, no reduction is known in any direction. Indeed, no such reduction is known 
even between the linear special cases, that is, generalized flows and minimum-cost circulations. In 
fact, the only known market setting captured by both is linear Fisher. Perfect price discrimination 
and ADNB are not known to be reducible to flows with convex objective. In contrast, spending 
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constraint utilities |40j are not known to be captured by our model, although they are captured by 
the other. 

As discussed in Section 14.41 it seems difficult to extend any generalized flow algorithm having 
separate cycle cancelling and flow transportation subroutines. While this includes the majority 
of combinatorial algorithms, there are some exceptions. Goldberg, Plotkin and Tardos ^12j gave 
two different algorithms: besides Fat-Path, they also presented MCF, an algorithm that uses a 
minimum-cost circulation algorithm directly as a subroutine. Hence for the concave setting, it could 
be possible to develop a similar algorithm using a minimum concave cost circulation algorithm, for 
example [18] or |24j as a black box. 

Another approach that avoids scaling is [43J for minimum-cost generalized flows and |33j for 
generalized flows: these algorithms can be seen as extensions of the cycle cancelling method, ex- 
tending minimum mean cycles to GAP's in a certain sense. While it does not seem easy, it might 
be possible to develop such an algorithm for concave generalized flows as well. 

In defining an e-approximate solution for the sink version of concave generalized flows, we allow 
two types of errors, both for the objective and for feasibility. A natural question is if either of 
these could be avoided. While the value oracle model as we use it, seems to need feasibility error, 
it might be possible to avoid it using a stronger oracle model as in |24] . One might also require a 
feasible solution as part of the input, as a starting point to maintain feasibility (For example if all 
lower bounds and node demands are and rjj(O) = on all arcs ij, then / = is always feasible). 
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